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Ammonia  borane  (AB)  is  a  candidate  material  for  on-board  hydrogen  storage,  and  hydrolysis  is  one  of  the 
potential  processes  by  which  the  hydrogen  may  be  released.  This  paper  presents  hydrogen  generation 
measurements  from  the  hydrolysis  of  dilute  AB  aqueous  solutions  catalyzed  by  ruthenium  supported 
on  carbon.  Reaction  kinetics  necessary  for  the  design  of  hydrolysis  reactors  were  derived  from  the  mea¬ 
surements.  The  hydrolysis  had  reaction  orders  greater  than  zero  but  less  than  unity  in  the  temperature 
range  from  16  °C  to  55  °C.  A  Langmuir-Hinshelwood  kinetic  model  was  adopted  to  interpret  the  data  with 
parameters  determined  by  a  non-linear  conjugate-gradient  minimization  algorithm.  The  ruthenium- 
catalyzed  AB  hydrolysis  was  found  to  have  activation  energy  of  76  ±  0.1  kj  mol“^  and  adsorption  energy 
of  -42.3  ±  0.33  kJ  mol“^  The  observed  hydrogen  release  rates  were  843mlH2min“^  (gcatalyst)“^  and 
8327  ml  H2  min“^  (gcatalyst)“^  at  25  °C  and  55  °C,  respectively.  The  hydrogen  release  from  AB  catalyzed 
by  ruthenium  supported  on  carbon  is  significantly  faster  than  that  catalyzed  by  cobalt  supported  on  alu¬ 
mina.  Finally,  the  kinetic  rate  of  hydrogen  release  by  AB  hydrolysis  is  much  faster  than  that  of  hydrogen 
release  by  base-stabilized  sodium  borohydride  hydrolysis. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Energy  issues  and  environmental  concerns  have  generated  inter¬ 
est  in  hydrogen  as  a  transportation  fuel.  On-board  hydrogen  storage 
poses  a  challenge  to  the  development  of  fuel  cell  powered  vehicles. 
Chemical  hydrides,  such  as  ammonia  borane  (AB),  are  a  plausible 
option  for  on-board  hydrogen  generation.  AB  is  a  non-toxic  material 
that  can  be  handled  safely  [1].  It  has  the  highest  material  hydro¬ 
gen  content  (about  19.6  wt.%)  among  all  amine  boranes.  AB  has  an 
energy  density  of  2.74  kWh  l“^  which  is  higher  than  the  energy 
density  of  2.36  kWh  1“^  of  liquid  hydrogen  [2].  However,  tempera¬ 
tures  above  500  °C  are  required  for  the  complete  dehydrogenation 
ofAB[2]. 

Dehydrogenation  through  AB  thermolysis  in  its  neat  form  [3-5] 
and  deposited  in  mesoporous  silica  [6]  have  been  reported  in  the 
literature.  Nucleation  seeding  of  AB  through  heat  treatment  of  the 
material  accelerates  hydrogen  release  rate  but  destabilizes  the  solid 
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fuel  [7].  The  effect  of  exothermic  hydrogen  release  on  stability  and 
safety  of  solid  AB  has  been  reported  in  recent  US  DOE  Laboratory 
studies  [7,8].  These  reports  emphasize  that  purity  is  critical  to  ther¬ 
mal  stability  of  solid  AB  [7,8].  It  is  reported  [8]  that  1  kg  of  AB, 
corresponding  to  1.91  of  closely  packed-pellets  (30%  void),  would 
be  required  to  meet  the  hydrogen  flowrate  requirements  at  peak 
vehicle  load.  Effective  methods  for  regeneration  of  AB  from  the 
dehydrogenation  products  are  being  sought  [4,9,10]. 

Aqueous  AB  is  stable  at  room  temperature  [3,11,12]  and  its  con¬ 
tact  with  catalysts  can  liberate  hydrogen  on  demand.  Catalytic  AB 
hydrolysis  generates  three  moles  of  hydrogen  from  one  mole  of  AB 
consumed  and  the  reaction  can  be  represented  as  follows: 

NH3BH3  +  2H20“-^''nH+  +  BO2  +  3H2  (1) 

This  is  an  exothermic  reaction  with  a  heat  of  reaction  of 
-156kJmol“^  and  the  products  of  this  reaction,  borate  ion  or 
boric  acid,  are  environmentally  benign  [13].  Hydrogen  release 
from  AB  hydrolysis  in  the  presence  of  transition-metal  catalysts 
(at  0.33  wt.%  AB),  non-noble  metal  catalysts  (at  1  wt.%  AB),  and 
solid  acids  (at  0.33  wt.%  AB)  has  been  investigated  [11-13]  at  room 
temperature  in  pursuit  of  a  low-cost  and  efficient  catalyst  for 
the  reaction.  However,  the  kinetic  parameters  of  these  reactions 
have  not  been  determined.  Among  the  non-noble  metal  catalysts, 
10  wt.%  cobalt  supported  on  alumina  (C0/AI2O3 )  was  found  to  be  the 
most  active  catalyst  [12]  and  among  the  transition-metal  catalysts. 
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Nomenclature 

A  pre-exponential  factor  (mol  (kg  catalyst)”^ ) 

Cab  concentration  of  NH3  BH  ( kmol  ) 

Cab.o  initial  concentration  of  NH3BH3  (kmol  m“^ ) 
dp  catalyst  particle  diameter  (p.m) 

£a  activation  energy  (kj  mol“^ ) 

^^ads  heat  of  adsorption  (kJ  mol“^ ) 

K  Langmuir  adsorption  equilibrium  constant 

(m^  kmol“^ ) 

Kq  Langmuir  adsorption  equilibrium  constant  at  a  ref¬ 
erence  temperature  (m^  kmol“^ ) 
k  reaction  rate  coefficient  (kmol  s“^  (kg catalyst)”^ ) 

kn  nth-order  reaction  rate  coefficient  (kmoF“^  s“^ 

(kg  catalyst)”^) 

meat  mass  of  catalyst  (mg) 

n  reaction  order 

Hab.o  moles  of  ammonia  borane  injected  into  the  batch 
reactor  (mol) 

nH2  moles  of  collected  H2  (mol) 

R  universal  gas  constant  (8.314J  mol“^  1<“^ ) 

A5°  entropy  change  of  the  Langmuir  adsorption 

isotherm  (J  mol“^  1<“^ ) 

T  reaction  temperature  (°C  or  K) 

To  reference  temperature  (K) 

Tw  wall  temperature  of  burette  (K) 

t  time  (s) 

Vsoi  volume  of  the  solution  (m^ ) 


20wt.%  platinum-on-carbon  (Pt/C)  produced  the  fastest  reaction 
rates.  It  has  been  observed  that  the  Pt/C  catalysts  are  highly  dis¬ 
persed  and  show  higher  catalytic  activity  than  Pt-black  alone  [11]. 
The  effect  of  Pt  catalyst  on  AB  hydrolysis  was  also  investigated 
in  Ref  [2].  The  activation  energy  (£a)  of  a  l<2Cl6Pt-catalyzed  AB 
hydrolysis  reaction  (at  0.13  wt.%  AB  in  D2O)  was  determined  to  be 
86.6  kj  mol“^  [2].  These  experiments  involved  a  temperature  range 
of  25-35  °C.  The  kinetic  parameters  reported  in  Ref.  [2]  involve  a 
second  order  dependence  of  the  AB  hydrolysis  rate  on  the  catalyst 
concentration  in  solution.  However,  a  significantly  lower  activation 
energy  =  62  kj  mol“^  for  cobalt-catalyzed  AB  hydrolysis  (1  wt.%) 
has  been  reported  for  a  temperature  range  of  20-40  °C  in  Ref. 
[12]. 

Recently,  a  preparative-scale  synthesis  of  AB  was  demonstrated 
in  Ref.  [1].  The  authors  examined  the  efficiency  of  a  series  of 
transition-metal  catalysts  such  as  PdCh,  NiCh.  C0CI2.  and  RUCI3; 
and  found  RUCI3  to  be  the  most  active  for  AB  hydrolysis.  They  also 
reported  liberation  of  ammonia  from  the  AB  hydrolysis  reaction 
at  relatively  high  AB  concentrations  (15  wt.%  and  25  wt.%).  Ammo¬ 
nia  liberation  was  not  detected  at  a  lower  (6  wt.%)  concentration 
[1]. 

A  complete  examination  of  the  effects  of  catalyst  choice  on  the 
reaction  rates  and  the  appropriate  reaction  models  and  parameters 
was  not  found  in  the  literature.  Based  on  this,  the  specific  objec¬ 
tives  of  the  present  study  are  to:  (1)  measure  the  intrinsic  kinetics 
of  Ru-catalyzed  AB  hydrolysis  under  isothermal  conditions,  after 
having  removed  the  effects  of  internal  and  external  diffusion,  using 
commercially  available  Ru  on  carbon  solid  catalyst;  the  kinetics 
model  developed  in  this  research  is  also  desired  to  interpret  the 
kinetics  measurements  obtained  in  concentrated  solutions  [14]; 
(2)  compare  the  effects  of  different  catalysts  on  the  kinetics  and 
discuss  these  effects  on  reactor  design  and  (3)  compare  the  AB 
hydrolysis  kinetics  with  that  of  sodium  borohydride  (SBH)  hydrol¬ 
ysis. 


The  experimental  methods,  the  kinetic  models,  the  method  for 
determination  of  the  model  parameters  and  conclusions  are  pre¬ 
sented  below. 

2.  Experimental  methods 

2.1.  Experimental  materials  and  conditions 

Ammonia  borane  (purity  >98%),  synthesized  by  a  Purdue  team 
[1  ],  and  commercial  3-wt.%  Ru  supported  on  cylindrical  carbon  pel¬ 
lets  (2  mm  in  diameter,  3  mm  in  length,  specific  area  ~1000  m^  , 
Alfa  Aesar),  ground  and  sized  prior  to  use,  were  used  in  the  AB 
hydrolysis  measurements.  In  past  hydrolysis  studies  in  this  labora¬ 
tory  [15],  this  form  of  supported  Ru  catalyst  was  found  to  be  more 
efficient  than  other  forms,  such  as  Ru  on  alumina  pellets,  and  Ru 
on  carbon  granules.  A  low  concentration  aqueous  AB  solution  (at 
1  wt.%)  and  a  small  amount  of  catalyst  (15.2  ±  0.1  mg)  were  used  to 
limit  the  reaction  rate  to  a  measurable  range  and  minimize  the  heat 
generation  during  the  experiments. 

Experiments  were  conducted  isothermally  at  four  tempera¬ 
tures:  16  °C,  25  °C,  35  °C  and  55  °C.  Catalyst  particle  size  of  22.5  pum 
and  stirring  speed  of  800  rpm  were  selected  to  eliminate  the  effects 
of  internal  and  external  mass  diffusion.  The  selections  of  catalyst 
size  and  stirring  speed  are  detailed  in  Section  2.4. 

2.2.  Experimental  apparatus 

The  experimental  apparatus  is  shown  schematically  in  Fig.  1.  A 
25  ml,  three-necked  reaction  flask  was  preloaded  with  the  ground 
catalyst  and  a  magnetic  stirring  bead.  The  ground  catalyst  was  pre¬ 
wetted  by  2  ml  of  deionized  water  to  eliminate  the  influence  of 
initial  wetting  [15].  The  reaction  flask  was  submerged  in  a  water 
bath  and  preheated  to  a  desired  reaction  temperature  using  a  ther¬ 
mostatic  circulator.  AB  (1  wt.%)  solution  was  prepared  in  another 
flask  and  preheated  to  the  desired  reaction  temperature  by  sub¬ 
merging  into  another  water  bath.  Preheating  the  AB  solution  is 
desired  to  minimize  the  heating  time  to  achieve  isothermal  con¬ 
dition  [15].  Approximately  2.53  millimoles  (around  7.8  ml)  of  the 
preheated  AB  solution  were  then  drawn  and  injected  into  the  reac¬ 
tion  flask. 

The  AB  solution  and  the  ground  catalyst  were  well  mixed 
by  using  the  magnetic  stirrer.  A  hypodermic  type-T  (copper- 
constantan)  thermocouple  with  a  stainless  steel  sheath  was  used  to 
record  the  temperature  of  the  solution  in  reaction  flask.  The  reac¬ 
tion  temperature  was  controlled  within  a  variation  of  ±0.6  °C  by 
the  thermostatic  circulator.  The  evolving  hydrogen  was  collected  in 
a  gas  burette  with  a  resolution  of  1.0  ml.  The  amount  of  collected 
hydrogen  in  the  burette  was  calculated  using  the  ideal  gas  law  with 
the  consideration  of  pressure  variation  caused  by  the  change  of 
water  column  height.  Since  the  mass  of  hydrogen  gas  evolved  from 
the  reactor  was  much  smaller  than  the  mass  of  the  burette  includ¬ 
ing  the  water  contained,  an  energy  balance  analysis  revealed  that 
the  differences  between  the  temperatures  of  the  collected  hydro¬ 
gen,  the  water  and  the  burette  were  less  than  0.1  °C.  As  a  result,  the 
burette  outer  wall  temperature,  Tw.  was  measured  to  represent  the 
accumulated  hydrogen  gas  temperature. 

2.3.  Experimental  uncertainties 

The  uncertainties  in  the  gas  temperature  and  gas  volume  mea¬ 
surements  caused  ±0.7%  experimental  uncertainties  in  the  mole 
of  hydrogen  collected  or  the  mole  of  AB  consumed.  Also,  ran¬ 
dom  errors  in  time,  meat  and  V^oi  measurements  caused  ±1.83% 
experimental  uncertainty  (with  95%  confldence)  based  upon  four 
repeated  tests  under  same  conditions.  As  a  result,  the  overall  exper¬ 
imental  uncertainties  in  the  mole  of  hydrogen  collected  or  the  mole 
of  AB  consumed  were  estimated  to  be  ±2.0%  with  95%  confldence. 
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Fig.  1.  Experimental  apparatus. 


2.4.  Selections  of  catalyst  particle  size  and  stirring  speed 

To  eliminate  the  internal  mass  diffusion  effects,  AB  hydroly¬ 
sis  using  different  catalyst  particle  sizes  was  studied  at  the  four 
experimental  temperatures.  The  catalyst  particles  were  ground  and 
sieved  into  different  size  ranges.  The  sieve  numbers  and  the  aver¬ 
age  catalyst  particle  sizes  (dp)  trapped  between  successive  sieves 
are  presented  in  Table  1.  At  each  temperature,  the  catalyst  particle 
size  for  which  the  measured  reaction  rates  were  independent  of 
the  particle  size  was  found.  The  lowest  particle  size  of  22.5  p.m,  at 
which  the  reaction  rates  were  found  to  be  free  from  internal  mass 
diffusion  effects  at  all  four  experimental  temperatures,  was  used  to 
obtain  kinetic  data. 

In  order  to  determine  a  suitable  stirring  speed  where  the  hydrol¬ 
ysis  reaction  would  be  free  from  external  mass  diffusion,  tests  were 
conducted  at  25  °C  and  55  °C  with  stirring  speeds  in  a  range  of  0  rpm 
to  1000  rpm.  The  fastest  reaction  rates  were  observed  at  a  stirring 
speed  of  800  rpm.  The  experiments  at  25  °C  are  shown  in  Fig.  2. 
The  rates  observed  at  1000  rpm  were  lower  than  those  at  800  rpm. 
We  assume  that  the  mass  transfer  rate  to  the  surface  of  the  catalyst 
decreases  as  a  result  of  reduced  relative  velocity  between  the  liquid 
and  the  catalyst  particles  caused  by  entrainment  at  the  highest  stir¬ 
rer  speed.  The  optimum  speed  of  800  rpm  was  selected  to  obtain 
kinetic  data  that  are  free  from  external  mass  transfer  effects. 

3.  Kinetic  models 

3.1.  nth-order  kinetics 

For  a  batch  reactor  with  a  volume  Vso\  and  a  catalyst  mass  meat* 
the  reaction  order  n  (n  ^  1 )  can  be  determined  using  the  following 
equation  ifa  plot  of  (C^gQ  -  C^g")/(1  -  n)  as  a  function  of  time  gives 
a  straight  line  through  the  origin  [15] 

rl-n  _  /-1-n  , 

'^AB,0  '^AB  _  {  rncaV<n\  ^ 

1-n  Vsol  )  ^  ^ 


where  Cab.o  is  the  AB  concentration  before  hydrolysis  and  Cab  is  the 
instantaneous  AB  concentration  during  the  reaction. 

3.2.  Langmuir-Hinshelwood  model 

For  a  liquid  phase  reaction  on  a  catalyst  surface  with  a  reac¬ 
tion  order  between  zero  and  one,  a  Langmuir-FIinshelwood  (LH) 
model  can  be  adopted  to  account  for  the  two  important  steps:  (a) 
equilibrated  adsorption  of  AB  molecules  on  the  surface  of  the  cat¬ 
alyst,  and  (b)  reaction  of  the  adsorbed  species  to  form  hydrogen. 
Following  Zhang  et  al.  [15],  the  reaction  rate  can  be  expressed  as 

d^AB  _  ^cat,  KCab  /on 

dt  ~  Vsoi  I+ZCCab  ^  ^ 

where  k  is  the  reaction  rate  coefficient  and  K  is  the  adsorption  equi¬ 
librium  constant.  The  reaction  rate  coefficient  can  be  expressed 
as 

k  =  >lexp(-^)  (4) 

And  the  adsorption  equilibrium  constant  can  be  expressed  as 

<=> 

where  Tq  =  298  K  is  the  reference  temperature. 

Eq.  (3)  represents  a  zero  order  reaction  when  KCab  is  much 
greater  than  unity  and  represents  a  first  order  reaction  when  /CCab 


Table  1 

Mesh  ranges  and  pore  sizes  of  the  sieves  used. 
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Fig.  2.  Effect  of  stirring  speed  on  hydrolysis  rate  at  25  °C,  the  time  axis  is  scaled,  as 
r  =  (('^cat/m*3j(n*g  Q/nAB,o))t,  to  account  for  differences  in  the  catalyst  masses  and 

moles  of  AB  injected  in  each  of  the  experiments.  q,  are  chosen  reference 

values  of  Hab.o  and  meat  in  any  one  of  the  experiments. 


Fig.  4.  Effect  of  solution  volume  on  hydrolysis  rate. 


is  much  less  than  unity.  Eq.  (3)  also  indicates  that  the  reaction  rate 
(dCAfi/df)  is  proportional  to  meat  and  inversely  proportional  to  Vgoi 
at  the  initial  part  of  the  reaction  when  /CCab  is  much  greater  than 
unity.  These  relationships  were  verified  with  a  range  of  meat  and 
Vsoi  in  preliminary  experiments.  With  all  other  reaction  conditions 
same,  AB  was  hydrolyzed  at  25  °C  with  meat  of  7.4  mg  (A),  15.2  mg 
(B)  and  29.4  mg  (C)  or  with  the  ratio  of  1.0:2.1:4.0.  Fig.  3  presents 
the  change  of  ( Cab, o -Cab)  with  reaction  time,  the  observed  rates  of 
reactions  A,  B,  and  C  had  a  ratio  of  1:2.1 :3.8.  Also,  with  all  other 
reaction  conditions  same,  AB  was  hydrolyzed  with  of  3.92  ml 
(A),  7.79  ml  (B)  and  12.75  ml  (C)  or  with  the  ratio  of  1.0:2.0:3.2.  As 
illustrated  in  Fig.  4,  the  observed  rates  of  reactions  A,  B,  and  C  had 
a  ratio  of  l:(l/2.0):(l/2.9). 

4.  Determination  of  Langmuir-Hinshelwood  kinetics 
model  parameters 

A  non-linear  approach  was  adopted  to  determine  the  four  kinet¬ 
ics  parameters.  A,  Fa*  f^o*  ^nd  AHads-  With  the  assumption  of  zero 
order  reaction  (infinite  K)  in  Eq.  (3),  (Cab.o  -  Cab)  is  a  linear  func¬ 


tion  of  t  and  the  slopes  of  these  straight  lines  at  four  experimental 
temperatures  provided  the  initial  guess  values  of  k.  Then,  an  Euler 
predictor-corrector  finite  difference  scheme  based  on  Eq.  (3)  was 
used  to  predict  the  change  of  AB  concentrations  at  various  time 
steps  at  four  temperatures  as  follows: 


^predicted  ^  ^cat , 

S+1  - 


KQ 


Vsol  1+/^Q 


At 


^TOrrected  _  q.  _ 


KC 


predicted 

i+l 


2V’sol  \  1  +  1  +  KCf 


predicted 

i+l 


At 


(6) 


Error  minimization  algorithms  were  then  used  to  estimate 
Kq  and  AHads  simultaneously  by  comparing  the  measured 
and  predicted  change  of  AB  concentrations.  Following  Pow¬ 
ell’s  conjugate-gradient  algorithm  in  conjunction  with  Brent’s 
line  minimization  algorithm  [16],  Kq  and  AH^ds  were  estimated 
simultaneously  as  a  vector  Y  with  initial  guesses,  for  example 
(2.5  m^  kmol“^ ,  -44  kj  mol“^ ),  and  two  minimization  directions  Hi 
(0.1  m^  kmol“^ ,  0)  and  U2  (0, 1  kJ  mol“^ ).  After  updating  K  with  the 
estimated  Kq  and  AHads*  the  k  values  at  the  four  different  tempera¬ 
tures  were  further  refined  by  minimizing  the  differences  between 


Fig.  3.  Effect  of  catalyst  mass  on  hydrolysis  rate. 


t,  min 


Fig.  5.  Aqueous  AB  hydrolysis  at  four  temperatures.  Reaction  conditions: 
Vsol  =  9-7  ml,  Hab.o  =  2.53  mmol,  meat  =  15.2  mg. 
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Fig.  6.  Reaction  orders  (n)  of  Ru-catalyzed  AB  hydrolysis  at  four  temperatures. 


Fig.  8.  Comparisons  of  measured  and  predicted  AB  concentration  changes  at  four 
temperatures. 


the  predictions  and  measurements.  Finally,  A  and  Fa  were  obtained 
from  the  Arrhenius  plot  using  updated  k  values.  This  procedure  was 
repeated  until  the  change  in  Kq  was  less  than  0.1  m^  kmol“^  and  the 
change  in  AHads  was  less  than  0.5  kj  mol“^  It  was  verified  that  the 
final  results  were  independent  of  the  initial  guess  values.  The  uncer¬ 
tainties  (bounded  limits)  in  the  kinetics  parameters  were  estimated 
by  varying  the  measured  Cab  data  by  ±2.0%. 

5.  Results  and  discussion 

The  measured  moles  of  hydrogen  evolved  per  mole  of  injected 
AB  (nH2/nAB,o)  ^re  presented  in  Fig.  5.  As  expected,  three  moles  of 
hydrogen  were  generated  for  one  mole  of  AB  consumed  at  the  end  of 
the  reaction.  Also,  the  reaction  rate  increases  rapidly  with  increas¬ 
ing  temperature.  Around  25  °C,  the  hydrogen  release  rate  from  AB 
hydrolysis  using  3  wt.%  Ru-C  is  843  ml  FI2  min“^  (g  catalyst)”^ .  This 
rate  is  13  x  faster  than  that  of  AB  hydrolysis  using  10  wt.%  C0-AI2O3, 
which  is  63mlFl2min“^  (gcatalyst)”^  [12].  Although  the  cost  of 
cobalt  is  much  less  than  that  of  ruthenium,  an  on-board  AB  hydrol¬ 
ysis  reactor  using  cobalt  will  be  much  bigger  and  therefore  heavier 
than  the  one  using  ruthenium  to  provide  the  same  desired  hydrogen 


flow  rate.  Considering  the  importance  of  reducing  weight  and  vol¬ 
ume  of  an  on-board  hydrogen  storage  system,  a  tradeoff  between 
the  reactor  size  and  cost  should  be  considered  and  the  final  choice 
of  the  catalyst  should  involve  an  optimization  of  these  two  aspects. 

It  is  also  found  that  the  Ru-catalyzed  AB  hydrolysis  reaction  had 
an  average  order  of  0.45  in  AB  in  the  temperature  range  of  16-55  °C. 
Different  guess  values  of  n  were  tried  until  linear  relationships 
between  (Cj^Q  -  Cj“")/(1  -  n)  and  t  were  achieved  as  shown  in 
Fig.  6. 

The  estimated  parameters  for  the  LH  kinetics  model  are: 
/Co  =  30.4  ±  5.6  m^  kmol"^ ,  AH^ds  =  -42.3  ±  0.33  kJ  mol"^ ,  A  = 
5  ±0.25  X  10^2  mol  s“^  (kgcatalyst)“^  and  £a  =  76±0.1  kJmol“^ 
The  Arrhenius  plot,  that  was  used  to  determine  A  and  Fa  using  the 
optimized  k  values  at  the  four  experimental  temperatures,  is  shown 
in  Fig.  7.  The  predicted  changes  of  AB  concentrations  with  time  at 
the  four  temperatures,  using  Eq.  (3)  and  the  estimated  parameters, 
are  depicted  in  Fig.  8.  Excellent  matches  were  achieved  at  16  °C  and 
35  °C  while  slight  differences  exist  at  25  °C  and  55  °C.  The  overall 
root  mean  square  error  is  3.9%  of  the  initial  concentration  of  the 
AB  solution. 

We  note  that  the  average  order  of  0.45  used  in  Fig.  6  is  con¬ 
sistent  with  the  LH  parameters  but  overlooks  the  change  in  order 


Fig.  7.  Arrhenius  plot  using  converged  reaction  rate  coefficients  at  four  tempera¬ 
tures. 


Fig.  9.  Changes  of  reaction  rate  coefficients  and  adsorption  equilibrium  constants 
with  temperature. 
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Table  2 

Reaction  properties  of  NH3BH3  and  NaBH4. 


AB 

SBH  [15] 

Fa  (kj  mol“U 

76 

66.9 

A(mols“^  (kgcatalyst)“U 

5  X  10^2 

2.19  X  lO^o 

kat25°C(mols“^  (kgcatalyst)“U 

0.25 

0.042 

AHads  (kJ  mol-i ) 

-42.3 

-35 

(Jmol-U<-U 

-141.8 

-130 

K at  25  °C  (m^  krnoP^ ) 

30.4 

220 

expected  over  the  wide  range  of  composition  at  each  temperature 
and  from  one  temperature  to  another.  We  regard  to  LH  model  as 
a  more  complete  description  of  the  kinetic  behavior  of  this  sys¬ 
tem. 

The  reaction  rate  coefficients  and  adsorption  equilibrium  con¬ 
stants  at  different  temperatures  are  depicted  in  Fig.  9.  The 
adsorption  equilibrium  constant  decreases  with  the  increasing 
temperature  moderately  indicating  stronger  adsorption  at  low  tem¬ 
peratures.  The  reaction  rate  coefficient,  however,  increases  with 
the  temperature  rapidly.  As  a  result,  the  reaction  rate  coefficient 
has  a  much  stronger  influence  on  the  overall  reaction  rate  than  the 
adsorption  equilibrium  constant. 

Also,  the  change  in  entropy  at  25  °C  of  the  adsorption  phe¬ 
nomenon  was  calculated  as  follows: 

A.0  AHads  -42,  300 ±330  (Jmor’) 

To  298.15  {K) 

= -141.8  ±1.1  Jmor’K-’  (7) 

The  negative  value  of  entropy  change  reflects  the  fact  that 
the  degree  of  randomness  decreases  in  an  adsorption  process,  as 
expected. 

Using  the  same  3  wt.%  Ru-C,  SBH  hydrolysis  kinetics  was  pre¬ 
viously  investigated  [15].  The  aqueous  SBH  solution  was  stabilized 
by  adding  NaOH.  Base-stabilized  aqueous  SBH  solution  had  a  pH 
value  of  14  which  is  higher  than  that  of  the  aqueous  AB  solu¬ 
tion  (9.1  [13]).  At  25 °C,  there  was  an  induction  time  observed  in 
SBH  hydrolysis  [15]  and  the  observed  hydrogen  evolution  rates  are 
134mlH2min“^  (g  catalyst)”^  and  843mlH2min“^  (g  catalyst)”^ 
for  SBH  and  AB  hydrolysis,  respectively.  At  55  °C,  the  observed 
hydrogen  release  rates  are  2047mlH2  min“^  (g  catalyst)”^  and 
8327  mlH2  min“^  (g  catalyst)”^  for  SBH  and  AB  hydrolysis,  respec¬ 
tively.  AB  hydrolysis  has  much  faster  kinetics  than  SBH  hydrolysis; 
so  a  packed-bed  AB  hydrolysis  reactor  would  be  much  smaller  and 
lighter  than  a  SBH  hydrolysis  rector  to  provide  identical  hydrogen 
flow  rates. 

The  reaction  properties  of  AB  and  SBH  hydrolysis  are  compared 
in  Table  2.  The  higher  value  of  /C,  at  25  °C,  implies  a  lower  reac¬ 
tion  order  in  case  of  SBH  hydrolysis  than  that  of  AB  hydrolysis  at 
room  temperature.  Also,  higher  Fa  and  AHads  of  AB  hydrolysis  indi¬ 
cate  stronger  temperature  dependences  in  both  reaction  (/<)  and 
adsorption  (/C)  processes  than  for  SBH  hydrolysis. 


6.  Conclusions 

The  following  conclusions  are  drawn  from  the  present  work: 

1.  A  non-linear  fitting  approach  based  on  Powell’s  conjugate- 
gradient  algorithm,  developed  in  this  study,  was  successful 
in  obtaining  reaction  rate  constants  with  the  Langmuir- 
Hinshelwood  (LH)  model. 

2.  Ruthenium  on  carbon  (Ru/C)  is  a  suitable  catalyst  for  the  ammo¬ 
nia  borane  hydrolysis  reaction.  After  normalizing  for  metallic 
catalyst  mass,  the  hydrogen  release  rate  of  AB  hydrolysis  is  13  x 
faster  when  using  Ru/C  catalyst  than  cobalt  catalyst  at  25  °C.  The 
effects  of  the  relative  rates  and  relative  costs  of  these  catalysts 
on  the  cost  of  ownership  need  to  be  investigated. 

3.  For  similar  reaction  conditions  and  a  lower  basicity  (pH),  hydro¬ 
gen  generation  from  AB  is  4x  and  6x  faster  than  hydrogen 
generation  from  sodium  borohydride  at  25  °C  and  55  °C,  respec¬ 
tively. 
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